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ABSTRACT 

We present new results on the nature of the Herbig-Haro-like object 73 (HHL 73, also known 
as [G84b] 1 1) based on narrow-band CCD Ha and [S II] images of the HHL 73 field, and Inte- 
gral Field Spectroscopy and radio continuum observations at 3.6 cm covering the emission of 
the HHL 73 object. The CCD images allow us to resolve the HHL 73 comet-shaped morphol- 
ogy into two components and a collimated emission feature of ~ 4-arcsec long, reminiscent 
of a microjet. The IFS spectra of HHL 73 showed emission lines characteristic of the spectra 
of Herbig-Haro objects. The kinematics derived for HHL 73 are complex. The profiles of the 
[S II] AA6717, 6731 A lines were well fitted with a model of three Gaussian velocity com- 
ponents peaking at Vlsr — —100, —20 and +35 km s -1 . We found differences among the 
spatial distribution of the kinematic components that are compatible with the emission from 
a bipolar outflow with two blueshifted (low- and high-velocity) components. Extended radio 
continuum emission at 3.6 cm was detected showing a distribution in close agreement with 
the HHL 73 redshifted gas. From the results discussed here, we propose HHL 73 to be a true 
HH object. IRAS 21432+4719, offset 30 arcsec northeast from the HHL 73 apex, is the most 
plausible candidate to be driving HHL 73, although the evidence is not conclusive. 

Key words: ISM: jets and outflows - ISM: individual: [G84b] 11, HHL 73, HH 379, 
IRAS 21432+4719 



1 INTRODUCTION 

Herbig-Haro-like objects (HHLs) appear as bright nebulosities in 
the red optical images. Spectroscopy of HHLs indicates that a frac- 
tion of them are reflection nebulae originated by the scattered light 
from a stellar source associated with the nebulosity, while only 
a fraction of HHLs are Herbig-Haro (HH) objects arising from 
shock-excit ed gas. HHL 73 , also known as [G84b] 11, was first 
reported by Gyulbudaghian (1984) as object 11 found in his sur- 
vey of DSS red plates. In these plates, the object shows a comet- 
shaped morphology, being located close to the western border of 
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a region with high visual extinction. In fact, HHL 73 is found as- 
sociated with a dark cloud complex in Cygnus, close to the Lynds 
Dark Cloud 1035A (LDN 1035A), and lying in a region with clear 
signs of recent star form ation activit y. The estimated distance to the 
region is 900 ± 100 pc dEliasl 19781) . 

The source IRAS 21432+4719 is located close to HHL 73. 
Both sources are engulfed within a clump of high-density molec- 
ular gas, traced by the emission of the NH3 and CS molecules, 
whose emissi on peaks around the HHL 73 and IRA S 21432+4719 
positions (see Anglada, Sepulveda & Gomez 1997 and references 
therein for a wider description on the high-density molecular gas 
emissions in the HHL 73 field). Emissio n from high-velocit y 
molecular gas in this field is also reported bv lDobashi et alj|l993l) . 
who detect a highly asymmetric CO bipolar outflow. The emis- 
sion of the redshifted outflow lobe is spatially more extended than 
the blueshifted lobe and encompasses the region around IRAS and 
HHL 73. In additio n, the detection of a water maser near HHL 73 
was r eported by (Gyulbudaghian, Rodrig uez & Mendoza-Torresl 
d 19871) . 

iDevine. Reipurth & Balivl i ll 9971) imaged the HHL 73 field 
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Figure 1. CCD H«-band image of the HHL 73 field obtained at the 
Isaac Newton Telescope (INT, ORM). Arrows point to the HHL 73 
and HH 379 objects. The position of IRAS 21432+4719 and the er- 
ror ellipse are marked. The tilted cross marks th e position given by 
iGvulbudaghian. Rodriguez & M endoza-To rresl i 19871) for the H2O maser. 
The image FOV is ~ 8.5 X 6 arcmin 2 . North is up and east is to the left. 



through narrow-band Ha and [S II] AA6717, 6731 A filters. They 
find two stellar-like condensations ~ 2.5 arcmin west of HHL 73 
with similar Ha and [S II] surface brightness, surrounded by dif- 
fuse emission. They designate their finding as the new Herbig-Haro 
object HH 379. In addition, they suggest that the exciting source 
of HH 379 may be embedded within the cometary nebula which 
appears in their images at a position close to IRAS 21432+4719 
(thus, at the position given bv lGvulbudaghianl Jl984l) for HHL 73). 
These authors also discuss a plausible association between the 
IRAS source and the undetected exciting source of H H 379. It 
should be noted that iDevine. Reipurfh & Ballvl d 19971) never re- 
fer tojhe_oj)tkalccmetary nebula a s the HHL 73 object re ported 
bv lGvulbudaghianl dl984» . Later on. lDavis et all d200lL[2003h per- 
formed near-//? echelle spectroscopy at the HHL 73 position. They 
detect emission from H2 A2.122 /im and [Fell] A 1.644 fim lines, 
while no continuum emission is detected at the HHL 73 position. 
Both lines show a rather complex profile, the width of the [ Fe II] 
line being roughly twice that of the H2 line. lDavis et al. ( 2003]) sug- 
gest that the [Fe II] emis sion may trace the higher- velocity compo- 
nent of ajet. In addition. iDavis et al. (2001, 2003 ) also suggest tha t 
the optical nebula found in the IDevine. Reipurth & Ballvl J 19971) 
images harbours the true exciting source of HH 379, thus discard- 
ing the possibiility that IRAS 21432+4719 was driving HH 379. 
This yet undetected exciting source of HH 379 is named as HH 
379-IRS by these authors. In order to help with the identification of 
the sources mentioned here, we show in Fig.Q]a CCD Ha image of 
the HHL 73 field where all these sources have been marked. 

With the aim of unravelling the nature of HHL 73, we included 
this target within a programme of Integral Field Spectroscopy (IFS) 
of Herbig-Haro objects using the Potsdam Multi- Aperture Spec- 
trophotometer (PMAS) in the wide-field IFU mode PPAK. Addi- 
tional CCD narrow-band images and radio continuum 3.6 cm data 
have been analysed to complement the study. We present in this pa- 
per the results on the kinematics and physical conditions of HHL 73 
obtained from this IFS pilot programme. 

The paper is organized as follows. The observations and data 
reduction of CCD imaging, radio continuum and IFS observations 
are described in § 2.1, 2.2 and 2.3 respectively. Results are given in 



§ 3: CCD imaging in § 3.1, radio continuum in § 3.2 and IFS in § 
3.3, including results from integrated emission (§ 3.4) and a line- 
profile decomposition model for [S II] profiles (§ 3.5). A global dis- 
cussion is made in 8 4. In 8 5 the main conclusions are summarized. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Narrow-band CCD images 

Narrow-band CCD images of the HHL 73 field were made on 2006 
August 5 at the 2.6 m Nordic Optical Telescope (NOT) (ORM, 
La Palma, Spain) using the Service Time facility. The images 
were obtained with the Andalucia Faint Object Spectrograph and 
Camera (ALFOSC). The image scale was 0.188 arcsec pixel - . 
The effective imaged field was ~ 5x5 arcmin 2 . Images were 
obtained through filters of Ha (central wavelength, A c = 6564 A, 
bandpass AA = 33 A, which mainly includes emission from this 
line, although some contamination from [Nil] emission is also 
present) and [S II] (A c = 6725 A, AA = 60 A), which includes 
the two red [S II] lines. The integration time was 1800 s in Ha and 
2400 s in [S II]. The seeing was 0.6-0.8 arcsec. In order to check 
whether there is significant continuum contribution to the HHL 73 
red emission, additional exposures (600 s integration time) were 
acquired through a narrow-band, off-line filter, but no continuum 
signal around the HHL 73 location was detected. All the images 
were processed with the standard tasks of the IRAlfl reduction 
package. Astrometric calibration of the two final images was 
performed. In order to register the images, we used the coordinates 
from the USNO Catalogu^l of eight field stars well distributed on 
the observed field. The rms of the transformation was 0.12 arcsec 
in both coordinates. 

In addition to the NOT images, we also show in Fig.fT]a CCD 
image of the HHL 73 field obtained on 2003 December 15 at the 
2.5 m Isaac Newton Telescope (INT) through a narrow-band filter 
of central wavelength A c = 6564 A and bandpass AA = 50 A that 
also includes emission from the [Nil] AA6548, 6583 A lines. The 
integration time was 3600 s. Astrometric calibration of this image 
was also performed with an accuracy of a fraction of the pixel size 
(0.24 arcsec rms in both coordinates). This image has a spatial res- 
olution significantly lower than that of the NOT images since the 
image scale was 0.54 arcsec pixel" 1 and, in addition, it was ob- 
tained under worse seeing conditions (^3 arcsec). However, the 
INT image is deeper than the NOT images and has been used in 
our study to look for faint emission nebulae in the field that might 
be related to the HHL 73 object and/or with HH 379. 



2.2 Radio continuum and H2O maser emission 

We observed the 3.6 cm continuum emission towards HHL 73 with 
the Very Large Array of the NRAC0 in the C configuration during 



1 IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 

2 The USNOFS Image and Catalogue Archive is operated by the United 
States Naval Observatory, Flagstaff Station. 

3 The National Radio Astronomy Observatory is operated by Associated 
Universities Inc., under cooperative agreement with the National Science 
Foundation. 
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Figure 2. Close-up of the HHL 73 images obtained at the Nordic Optical Telescope (NOT, ORM) through the line filters of Ha (left) and [S II] (right). In 
both images, the spatial structure of HHL 73 emission is resolved into two components. A microjet of ~ 4 arcsec is also visible in the two band images. The 
FOV of each image is ~ 50 X 42 arcsec 2 . North is up and east is to the left. 



1997 August 14. The phase centre was located at RA( J2000) = 
21 h 45 m 10:46; Dec.(J2000) = +47° 33' 21'.' 2. The phase cali- 
brator used was 2202+422, and the bootstrapped flux density ob- 
tained was 3.69 ± 0.07 Jy. The absolute flux density calibrator 
was 00137+331 (3C48), for which a flux density of 3.15 Jy was 
adopted. Clean maps were obtained using the task MAPPR of AIPS 
with the robust parameter set equal to 5, which is close to natural 
weighting. Since the signal-to-noise ratio of the longest baselines 
was low, we applied a [u, u)-taper function of 50 kA. The resulting 
synthesized beam size was 4.7 x 4.3 arcsec 2 (P. A. = —78°), and 
the achieved la noise level was 20 /iJy beam -1 . 

VLA archive observations of the water maser line at 
22.2351 GHz (6ie-523 transition) carried out in the B configuration 
during 1994 July 4 were analysed. The phase centre was located at 
RA(J2000) = 21 h 45 m 08 ! :37; Dec.(J2000) = +47°32'49'.'l. 
The bootstrapped flux obtained for the phase calibrator 2202+422 
was 3.08 ± 0.09 Jy, and the adopted flux density for the absolute 
flux calibrator (3C48) was 1.13 Jy at 1.3 cm. The observations were 
made using the 2IF mode, with a bandwidth of 6.30 MHz with 31 
channels of 203 kHz resolution (2.66 km s -1 at 1.3 cm) centred at 
—2 km s , plus a continuum channel that contains the central 75 
percent of the bandwidth. The synthesized beam was 0.36 x 0.31 
arcsec 2 (P. A. = —11°), and the lcr noise level per spectral channel 
achieved with natural weight was 8 mJy beam - . 



2.3 Integral Field Spectroscopy (IFS) 

Observations of the HHL 73 object were made on 2004 November 
22 with the 3.5 m telescope of the Calar Alto Observatory (CAHA). 
Data were acquired with the Integral Field Instrument Pots dam 
Multi-Aperture Spectrophotometer PMAS teothetalj|2005h us- 
ing the PPAK configuration that has 331 science fibres, covering 
a hexagonal FOV of 74x65 arcsec 2 with a spatial sampling of 2.7 
arcsec per fibre. In addition, PPAK comprises 36 fibres to sam- 
ple the sky, located ~ 90 arcsec from the centre, packed in small 
bundles of six fibres each one, and distributed at the edges of the 
science hexagon. Fifteen calibration fibres are interpolated in be- 
tween the science and sky ones in the pseudo-s lit, and connecte d 
directly to the PMAS internal calibration unit feelz et alj|2006h . 
The 11200 grating was used, giving an effective sampling of 0.3 A 
pix -1 (~ 15 km s -1 for Ha) and covering the wavelength range 



~ 6500-7000 A, thus including characteristic HH emission lines 
in this wavelength range (Ha, [N II] AA6548, 6584 A and [S II] 
AA6717, 6731 A). The spectral resolution (i. e. instrumental profile) 
is ~ 2 A FWHM (~ 90 km s -1 ) and the accuracy in the determina- 
tion of the position of the line centroid is ~ 0.2 A (~ 10 km s _1 for 
the strong observed emission lines). We made a single pointing of 
1800-s exposure time centred on the SIMBAD position of HHL 73 
(RA(J2000) = 21 h 45 m 08n, Dec.(J2000) = +47°33' 6") cov- 
ering the whole optical emission of the object (~ 30 arcsec). 

Data reductio n was performe d using a preliminary version of 
the R3D software iSanchezj|2006h, in combination with IRAF and 
the Euro3D packages (Sanchez 2004). The reduction consists of the 
standard steps for fibre-based integral field spectroscopy. A mas- 
ter bias frame was created by averaging all the bias frames ob- 
served during the night and subtracted from the science frames. 
The location of the spectra in the CCD was determined using a 
continuum-illuminated exposure taken before the science expo- 
sures. Each spectrum was extracted from the science frames by 
co-adding the flux within an aperture of five pixels along the cross- 
dispersion axis for each pixel in the d ispersion axis, and stored it 
in a row-stacked-spectrum (RSS) file (Sanchez 2004). At the date 
of the observations we lacked a proper calibration unit for PPAK. 
Therefore, the wavelength calibration was performed iteratively. 
First, we illuminated the science fibres using dome strahler lamps 
with the telescope at the location of the objects. These lamps have 
well defined emission lines in the wavelength range considered, al- 
though the line elements and wavelength are ignored at this step. 
These exposures were used to correct for distortions in the wave- 
length solution, fixing all the spectra to a common (ignored) so- 
lution. After that, the wavelength solution was found using HgNe 
lamps that illuminated only the calibration fibres. Differences in 
the fibre-to-fibre transmission throughput were corrected by com- 
paring the wavelength-calibrated RSS science frames with the cor- 
responding continuum illuminated ones. The final wavelength so- 
lution was set by using the sky emission lines found in the ob- 
served wavelength range. The achieved accuracy for the wave- 
length calibration was better than ~ 0.1 A (~ 5 km s _1 ). Ob- 
servations of a standard star were used to perform a relative flux 
calibration. A final datacube containing the 2D spatial plus the 
spectral information of HHL 73 was then created from the 3D 
data by using Euro3D tasks to interpolate the data spatially until 
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Figure 3. Left panel: Continu um emission at 3.6 cm (white contours) overlai d on the A'-band 2MASS image (grey scale). The cross marks the position of the 
H2O maser detected by Gyulbudaghian, Rodriguez & Mendoza-Torres 1 1987). The synthesized beam is shown at the bottom left corner, and the black ellipse 
indicates the IRAS error ellipse. Right panel: Detail of the 3.6 cm continuum emission towards HHL 73 (white contours) overlaid on the [S II] image (grey 
scale). In both panels, contours levels are —5, —3, 2, 3, 5, 7, and 9 times the rms of the map, 20 p,3y beam . 



reaching a final grid of 2.16 arcsec of spatial resolution and 0.3 
A of spectral resolution. Further manipulation of this datacube, de- 
voted to obtaining integrated emission line maps, channel maps and 
position-velocity maps, were made using several common-users 
tasks of STARLINK, IRAF and GILDAS astronomical packages an d 
IDA jGarcfa-Lorenzo. Acosta-Pulido & Me gias-Ferna ndej l2002h . 
a specific suit of IDL software for analysing 3D data. 



3 RESULTS 

3.1 Narrow-band CCD images 

Figure Q] shows the Ha image of the HHL 73 field obtained at the 
INT. The HHL 73 object appears projected on a dark cloud near 
its southwestern border. The object presents in this Ha image the 
comet-shaped morphology already found in the DSS plates. Be- 
cause of the low effective spatial resolution of the image result- 
ing from both the pixel scale and the seeing, we were unable to 
resolve any structure in the morphology of the Ha emission. Nei- 
ther were we able to elucidate from this image whether the emis- 
sion of HHL 73 originates mainly from scattered starlight from an 
embed ded point-source, as suggested bv lDevine, Reipurth & Bailvl 
(1997), or whether the emission mainly arises from gas photoion- 
ized or shock-excited by an external source. The error ellipse of 
IRAS 21432+4719 does not encompass HHL 73. From the astrome- 
try performed on the Ha image, we derived that the IRAS catalogue 
position is ~ 30 arcsec northeast of the HHL 73 Ha emission peak 
(~ 0.13 pc for a distance of 900 pc). A 2MASS source is found 
very close to HHL 73 (there is an offset of ~ 1.5 arcsec between 
the position given for the infrared source in the 2MASS catalogue 
and the position of the Ha emission peak of HH L 73). The neb- 
ulosities cited bv lDevine. Reipurth & Ballvl d 19971) as HH 379 are 
also visible in our INT image close to a bright field star (see Fig. 
[3. However, despite a careful inspection of the HHL 73 field im- 
aged in Ha, we were unable to identify any other nebular emission 
suggestive of being part of a putative optical jet related to HH 379. 

The narrow-band Ha and [S II] images obtained at the NOT 
have less depth than the INT image. However the higher spatial res- 
olution achieved in the NOT images make them suitable for resolv- 
ing the structure of the HHL 73 object. Fig. [2] shows the close-up 
of the NOT images including the HHL 73 object in the Ha (right) 



and [S II] (left) bands. As can be seen from the figure, these images 
allow us to resolve several components in the emission of HHL 73: 
the comet-shaped morphology detected in the optical INT images 
of lower spatial resolution now appears split into two components 
of different surface brightness with a darker lane between them. 
The estimated distance between the peaks of the two HHL 73 com- 
ponents is ~ 2.5 arcsec (~ 0.01 pc for a distance of 900 pc). In ad- 
dition, a lower brightness and more collimated emission tail, show- 
ing a couple of brightness enhancements suggestive of the "knots" 
of a microjet, emanates from the brighter northwestern component 
and extends ~ 4 arcsec towards the west at P. A. = 267°. All 
these HHL 73 components are found in both NOT images. How- 
ever, some differences are appreciated when the emission from the 
two bands is compared: the northwestern component is the brighter 
one in both filters, Ha and [S II]. The southeastern component ap- 
pears dimmer, as compared with the northwestern one, in Ha than 
in [S 11]. The emission tail is brighter and has higher collimation in 
[S II] than in Ha (as would be expected from a jet). Furthermore, 
while the emission peak positions of Ha and [S II] are coincident 
in the southeastern component, the Ha peak position in the north- 
western component is offset by ~ 0.8 arcsec towards the southwest 
from the [S II] emission peak position. Thus, differences in the ex- 
citation conditions between the two HHL 73 components might be 
present. 

Finally, note t hat the position angle of HH 379, as seen from 
HHL 73, is 262° jDavis et al.ll200lU2003l) . close to that of the 
HHL 73 microjet. Thus, although the deep INT Ha image does 
not show any other nebulosity connecting both emissions, HH 379 
could form part of the HHL 73 optical outflow, provided that 
the outflow direction undergoes a bending of ~ 5°, as seen in 
other optical outfl ows propagating in a inhomogeneous medium 
jLopezetalJ 1995b . 



3.2 Centimetre continuum and H2O maser emission 

In order to complement the CCD and IFS data of HHL 73, we anal- 
ysed the data of the emission at radio wavelengths available for the 
HHL 73 field. 
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Table 1. Parameters of the continuum sources detected in the HHL 73 region. 





Position 


Flux density 


Deconv. size 


PA. 


Source 


a(J2000) <5(J2000) 


(mjy) 


(arcsec) 


(deg) 


VLA 1 (HHL 73) 
VLA2 


21 h 45 m 08!15 47°33'05'.'6 
21 h 45 m 12!77 47°33'43('6 


0.26±0.02 
0.33±0.03 


3.6x2.0 
7.5x0.0 


170 

9 



a CoiTected for primary beam response. 

3.2.1 Centimetre continuum emission 

Figure [3] (left panel) shows the continuum map at 3.6 cm towards 
HHL 73 overlaid on the K-bwiA 2MASS image. We detected two 
centimetre sources, VLA 1 and VLA 2, separated ~ 47 arcsec 
(0.2 pc in projection at the distance of the source), located near the 
catalogue position of IRAS 21432+4719. Both centimetre sources 
are also associated with 2MASS point sources. In Table 1 we list 
the positions, flux densities and deconvolved sizes of the detected 
sources in the HHL 73 region. In Fig. [3] (right panel) we show the 
3.6 cm centimetre continuum emission arising from the HHL 73 
object overlaid on the [S II] image. The peak position of the cen- 
timetre source VLA 1 is close (~ 0.7 arcsec south) to the posi- 
tion of the [S II] emission peak. The centimetre emission from the 
HHL 73 object is slightly elongated towards the secondary peak 
detected in the [S Il]-band images to the southeast. 

3.2.2 H2O maser emission 

In 1985 H2O maser emission was detected towards HHL 73 by 
iGvulbudaghian. Rodriguez & Mendoza-Torresl dl987l) using the 37 
m radio telescope of the Haystack Observatory. The H2O maser 
was located 10.2 arcsec west of the HHL 73 object, with a posi- 
tional uncertainty that can be estimated at ~15 arcsec, making it 
difficult to ascertain whether or not the maser is associated with 
HHL 73. In the archive observations analysed here we did not de- 
tect H2O maser emission towards HHL 73 above a 6a limit, in- 
dicating that at the epoch of observation (1994 July) the maser 
was out. The upper limit estimated for the peak intensity of the 
H2O maser emission was 48 mjy beam -1 , corresponding to a main 
beam brightness temperature upper limit Tmb = 1000 K. 

3.3 IFS: Integrated line-emission maps 

3.3.1 Emission line images 

We detected emission from characteristic HH lines within the ob- 
served wavelength range (Ha, [Nil] AA6548, 6584 A and [S II] 
A6717, 6731 A). Narrow-band images for these lines were gener- 
ated from the datacube. For each position, the flux of the emission 
line was obtained by integrating the signal over the whole wave- 
length range covered by the line and subtracting a continuum ob- 
tained from the wavelengths adjacent to the emission line. 

The narrow-band maps of HHL 73 obtained are shown in Fig. 
[4] As can be seen from the panels of the figure, the emission has 
a similar morphology in all the mapped lines. However it should 
be noted that the Ha emission is slightly more extended (~ 3 arc- 
sec) towards the south than the emission from the other lines. Fur- 
thermore, the [N II] emission appears more compact than the Ha 
and [S II] emission. This is mostly because the [N II] emission has 
a lower surface brightness, and the collimated emission extending 
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Figure 4. Integrated maps of the HHL 73 emission obtained from the 
datacube by integrating the signal within the wavelength range including 
the corresponding line labelled in each panel (Ha: A6557. 7-6569.1 A, 
[Nil]: A6579.6-6585.9 A, [S ll]A6716: 6712.1-6719.3 A and [Sn]A6731: 
6728.0-6733.8 A). Fluxes are displayed on a logarithmic scale, in units of 
10~ 1B erg s _1 cm -2 . The offsets of the spatial scale are relative to the 
position of the emission peak. 



towards the west of the HHL 73 main body is hardly detected in the 
[N II] lines. 

It should also be noted that the overall morphology derived 
from the Ha and [S II] IFS images is in good agreement with the 
morphology found from the INT CCD image. In particular, the IFS 
images also show the low-brightness and more collimated emis- 
sion feature originating west of the HHL 73 main body. However, 
the lower spatial resolution of the IFS images did not allow us to 
resolve the two components of the HHL 73 main body discovered 
from the Ha and [S II] NOT images. 



3.3.2 Kinematics 

In addition to the integrated flux, we derived for each position the 
flux-weighted mean radial velocitjQ (V), and the flux-weighted 
rms width of the line (AV 2 ) 1 / 2 , proportional to the FWHM, which 
are given by 



4 All the velocities in the paper are referred to the local standard of rest 
(LS R) frame. A Vlsr. = +3.6 km s — 1 for the parent cloud has been taken 
from lDobashi et alj Jl993ri . 
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Figure 5. Flux-weighted mean velocity maps. Velocities (km s _1 ) are re- 
ferred to the Vlsr °f tne parent cloud. The offsets of the spatial scale 
are relative to the position of the emission peak. Contours of the inte- 
grated flux of the corresponding line have been overlaid. The lowest contour 
level is 7 x 10' " 
erg s _1 cm -2 . 



erg s 1 cm 1 and the contour spacing is 3 X 10 



(V) = j / IvVd\ 



and 



AV 2 



where 



Iv(V- (V)fdV, 



1= IvdV 



(1) 



(2) 



(3) 



In spite of the spatial extension of the emission being quite com- 
pact, the velocity field derived for (V) appears rather complex for 
all the lines and shows appreciable changes through HHL 73, as 
can be seen from Figure [5] 

For all the mapped lines, the gas appears blueshifted around 
the emission peak, with (V) values of —20 km s _1 for Ha and 
[Nil], and —2.0 km s _1 for the [S II] lines. However, the highest 
blueshifted velocities are found ~ 5 arcsec west of the emission 
peak position, where values of —50 km s _1 for Ha and [N II], and 
—32.0 km s _1 for the [S II] lines have been derived. 

The FWHM values obtained from the (AV 2 ) 1/2 , corrected 
from the instrumental profile, are high for all the mapped lines. 
Around the emission peak position, FWHMs are ~ 190 km s _1 
for the Ha line and ~ 100 km s' 1 for the [Nil] and [S II] lines. 
These values are representative of the FWHMs found in the whole 
emitting region, although the FWHMs slightly increase towards the 
southeast and decrease towards the west from the emission peak po- 
sition. Such large line widths might be indicating turbulent gas mo- 
tions, which are more important close to the position proposed by 
iDavis et all j200lLl20"03h for an unseen embedded exciting source. 
Another possibility is that the high FWHM values derived were 
due to the contribution of several kinematic components of the line 
profiles. We shall analyse this possibility later by using a multicom- 
ponent line profile model to fit the [S II] profiles. 
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Figure 6. Ha channel maps centred on the VLsr radial velocities (km s —1 ) 
labelled in each panel. The spatial scale is given in one of the panels. 
Fluxes have been displayed on a logarithmic scale, in units of 10~ 16 
erg s _1 cm -2 . (Channels corresponding to +21.7 ^ Vlsr. ^ +62.9 
are affected by residuals of an inaccurate subtraction of the sky background 
at the north of the emission). 



3.3.3 Physical conditions 

Line-ratio maps of the integrated line fluxes are suitable for getting 
spatial information about excitation conditions (from [S Il]/Ha, 
[Nll]/Ha) and electron density (from [S II] 6717/6731) of the 
emitting gas. For the entire region covered by the compact bright 
HHL 73 emission, the derived [S II] 6717/6731 line ratios go from 
~ 0.6 to ~ 0.9, indicating a range of electron densities n c = 900- 
4000 cm -3 , the highest values being reached around the emission 
peak position. Electron density lowers towards the west in the low- 
brightness collimated emission tail (e. g. [S II] 6717/6731 line ratio 
~ 1, corresponding to n c ~ 600 cm -3 , is derived at a position 
~ 5 arcsec west of the peak position). These values were derived 
for T c = 10 4 K. 

The derived [S Il]/Ha, [N Il]/Ha line ratio distributions show 
a similar pattern. Both line ratios peak around the emission peak 
position, with [Sll]/Ha ~ 0.9 and [Nll]/Ha ~ 0.4. Line ratios 
decrease steeply (excitation increases) southeast of the peak, while 
the decrease is slower to the west along the low-brightness emission 
tail. 



3.4 IFS: Channel maps 

An inspection of the HHL 73 spectra at different positions revealed 
that the line profiles, being broad and asymmetric, display a variety 
of features suggesting complex gas kinematics. 
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Figure 7. Same as Fig.|6] but for [N II] A6584 , 
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Figure 8. Same as Fig. [6] but for the [S II] 6717 A emission line. 



We will now analyse the behaviour of the emission as a func- 
tion of velocity. First, we obtained the Ha, [N II] and [S II] channel 
maps to explore the 3D structure of the physical conditions as a 
function of the radial velocity. Furthermore, to analyse the com- 
plex line profiles found along the HHL 73 emission, we extracted 
spectra for each spaxel within the area covered by the emission and 
we fitted the line profiles through a multicomponent line profile 
model. 



3.4.1 Spatial distribution of the emission 

A more detailed sampling of the gas kinematics was obtained by 
slicing the datacube into a set of velocity channels along each emis- 
sion line profile, thus providing a "tomograph" of the emission 
line, not restricted to the centroid of the line, but including also the 
wings. Each slice was made for a constant wavelength bin of 0.3 A, 



Q 




Figure 9. Position- Velocity diagram of the Ha emission for a cut along 
the Dec. axis at RA ~ 6.5 arcsec west of the position of the intensity peak. 
The AS = position is marked by the dashed line. The contour spacing is 
10 per cent of the peak intensity. 




V LSR ( km s ) 

Figure 10. Position- Velocity diagram of the Ha emission for a cut along 
the RA axis at Dec. ~ 2.2 arcsec south of the position of the intensity peak. 
The Aa = position is marked by the dashed line. The contour spacing is 
10 per cent of the peak intensity. 



corresponding to a velocity bin of ~ 14 km s~\ The channel maps 
obtained for Ha, [N II] A6584 A and [S II] A6717 A are shown in 
Figs |6]|7] and [8] respectively. As we used the same wavelength bin 
for all the lines, the central velocity of a given channel map differs 
slightly from a line to another. 

As can be seen from Figs [6] and [8] the spatial distribution of 
the Ha and [S II] emission found from the line-integrated IFS maps 
(i. e. the compact, bright region that extends towards the west in a 
low-brightness and collimated tail) is also found in all the veloc- 
ity channel maps ranging from blueshifted Vlsr velocities of —30 
km s _1 up to redshifted velocities of +20 km s _1 in both emission 
lines. 

The emission spreads over a wide velocity range for all the 
mapped lines, the span of the velocity range being slightly differ- 
ent for each line. The Ha emission extends over the widest range 
as compared with the other lines, covering a range ~ 400 km s _1 
wide, from —185 to +200 km s , The strongest and more ex- 
tended Ha emission is found for the velocity channels ranging from 
-50 to +20 km s~\ The [S II] emission is detected within a veloc- 
ity range ~ 250 km s _1 wide, from —145 to +100 km s~ . The 
strongest and more extended [S II] emission is found for the veloc- 
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ity channels ranging from —50 to +45 km s 1 . The [N II] emission 
is detected within a velocity range ~ 300 km s" 1 wide, from -170 
to +105 km s _1 . The [N II] emission presents a similar spatial ex- 
tension in most of the velocity channels since the low-brightness 
emission tail is hardly detected in the [N II] line. 

Furthermore, we performed position-velocity (PV) diagrams 
at different right ascension and declination positions, both in the 
Ha and [S II] lines. Two of these Ha cuts at selected positions are 
shown in order to illustrate more clearly a trend found in the veloc- 
ity distribution of the emission. 

Figure [9] shows the PV diagram obtained for a cut at RA ~ 
6.5 arcsec west of the position of the emission peak (marked in the 
figure by a dashed line). This cut along the declination axis has been 
done at an RA position that intersects the beginning of the colli- 
mated emission tail. The figure clearly shows that the overall emis- 
sion appears blueshifted north and south of the emission peak, at 
Vlsr — —110 km s~ . Another velocity component is also found 
peaking at ~ —70 km s~ . Furthermore, note that the spatial dis- 
tribution of the velocities is slightly asymmetric, with the range of 
velocities increasing from north to south: e. g. the velocities range 



from —130 to —15 km s 



at 



2 arcsec north of the emission 



peak, while the velocities range from —110 to —45 km s~ at ~ 2 
arcsec south of the emission peak. 

FigurefTOlshows the PV diagram obtained for a cut at Dec. ~ 
2.2 arcsec south of the position of the emission peak. For this 
cut along the RA axis, the figure shows that there is also a clear 
asymmetry in the spatial distribution of velocities: the emission at 
blueshifted velocities is spatially more extended towards the west, 
while at redshifted velocities it is is more extended towards the east 
of the peak emission position. Two blueshifted velocity compo- 
nents, peaking at ~ —70 and ~ —25 km s~\ are also detected 
for this cut. 

The inspection of different PV diagrams obtained at several 
positions up to ^ 10 arcsec around the emission peak revealed 
a general trend, both in the Ha and [S II] lines, in the highest 
blueshifted velocities appearing towards the west of the emission 
peak position, with the velocity values increasing from west to east 
along HHL 73. This signature of a two- velocity distribution in the 
emission can be hardly appreciated in the channel maps. 

3.4.2 Density and excitation conditions 

We also explored the behaviour of the density and excitation as a 
function of the velocity by generating the line ratio channel maps 
for a velocity range from —135 to +1 10 km s" 1 . 

Regarding the kinematics of the emission, the [S II] 6717/6731 
line-ratio channel maps indicate that n c decreases with increasing 
velocities: the blueshifted gas has higher n c than the redshifted gas. 
This trend is found at different position along the HHL 73 emission. 
Around the position of the emission peak, we derived n c ~ 2200 
cm -3 for the —80 km s _1 velocity channel, while n c ~ 740 cm -3 
is derived for the —55 km s _1 velocity channel. 

Regarding the spatial distribution of the emission, the [S II] 
6717/6731 line ratio channel maps indicate that n c decreases from 
east to west for all the velocity channels, the n c variations being 
different depending on the velocity. In fact, the variations found 
for n e along the east-west axis seem to decrease with increasing 
velocity: e. g. in the —80 km s _1 channel, n c changes from ~ 4200 
to ~ 740 cm" 3 (i. e. a factor of ~ 5) as we move ~ 10 arcsec from 
east to west; in contrast, n c only changes from ~ 1400 to ~ 490 
cm -3 for the same spatial displacement in the +82 km s _1 channel. 

The [S Il]/Ha line-ratio channel maps show an increasing 
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Figure 11. Spectrum extracted at the emission peak position. Flux has been 
normalized to the flux of the Ha peak. The characteristic HH emission lines 
(Ho, [N II] A6548, 6583 A and [S II] A6717, 6731 A) are clearly seen. 
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Figure 12. Two-dimensional distribution diagram of the observed spec- 
tra in the spectral window 6553.1-6590.7 A, including the Ha and [Nil] 
A6584 A lines. Each spectrum has been obtained by averaging the spectra 
within an aperture of 2 X 2 spaxels (~ 4 X 4 arcsec 2 ) and its intensity has 
been normalized to the Ha peak. Spectra are overlaid on the Ha integrated 
map obtained as in Fig. [4] 



trend (i.e. decreasing excitation) with increasing velocity, from 
blueshifted velocities up to a velocity of ~ +40 km s _1 . For 
higher redshifted velocities, line-ratio values decrease (i.e. exci- 
tation increases). Thus the behaviour derived from the line ratios is 
that the excitation is higher at high absolute velocity values (both 
blueshifted and redshifted) than at low (absolute) velocity values. 

Regarding the spatial distribution, the general trend found for 
each velocity channel is a decrease in the line ratios (i. e. an increase 
of excitation) from east to west. 

3.5 IFS: Line-profile multicomponent analysis 

Figure fTTI shows the spectrum extracted at the Ha emission peak 
position. As can be seen in the figure, the spectrum is closely remi- 
niscent of the characteristic shock-excited HH spectra. Note in ad- 
dition that the line profiles are broad and asymmetric, suggesting 
that multiple kinematic components are contributing to the line pro- 
file. Furthermore, the inspection of the spectra extracted at different 
positions over HHL 73 reveals that the line profiles display a variety 
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Figure 13. Same as Fig. 1131 but for the [S II] emission. The spectral window 
selected for this plot is 6706.4-6739.5 A. 

of features. Around the emission peak, line profiles are asymmetric 
and broad, with a faint blue double peak suggesting the presence 
of several gaseous systems. To the north and south of the emission 
peak, the profiles clearly have contributions from at least two com- 
ponents, showing double peaks and/or broad wings. Towards the 
west of the emission peak, in the low-brightness collimated emis- 
sion tail, the line profiles are more complex, showing blue and red 
shoulders as well as broad wings. All the emission lines in the ob- 
served spectral range show similar features to those described here, 
as can be seen in Figs[T2]and[T3] where the spectra for the spectral 
ranges which include the Ha and [N II] A6584 A (Fig.l 12b and the 
[S II] AA6717, 6731 A emission lines (Fig. 1 13ft have been overlaid 
on the flux maps obtained by integrating the signal within the same 
spectral range. 

The complex line-profile pattern appears in both the Ha and 
in the [S II] lines. We used the [S II] lines to perform a line-profile 
decomposition to reproduce the profiles of the HHL 73 spectra, 
since the [S II] emission should be less contaminated from the emis- 
sion of a putative, undetected source and thus, should better trace 
the kinematics of the shock-excited outflow gas. The [S II] line 
profiles were well reproduced by assuming the contribution from 
three kinematic components. The line-profile decomposition was 
obtained from a model including three Gaussian components. The 
Gaussian fitting was performed using the DIPSO package of the 
STARLINK astronomical software. The wavelength difference be- 
tween the two [S II] lines was fixed according to their correspond- 
ing atomic parameters. Furthermore, the same FWHM of the Gaus- 
sian in each of the two [S II] lines was assumed for each component 
of the profile. 

According to the fit parameters, we identified three Gaussians 
with three different kinematic components well-separated in veloc- 
ity (wavelength), as shown in Fig. [14] two of the velocity com- 
ponents are blueshifted and the third is redshifted. In the follow- 
ing, we refer to the blueshifted high-velocity component as BHVC 
(filled triangles in the figure), to the blueshifted low-velocity com- 
ponent as BLVC (open dots) and to the redshifted velocity com- 
ponent as RVC (filled dots). Reasonably good results were ob- 
tained from this three-Gaussian components model along the whole 
HHL 73 emission (see Fig. Q3] for some examples, where the fits 
obtained at selected positions and their residual are plotted). The 
Vlsr velocities of the BHVC range from —140 to —90 km s _1 , 



L_i i i I i i i I i i i I i_J 

6712 6714 6716 6718 

^central 

Figure 14. Central wavelength vs. FWHM for all the lines obtained by 
fitting the HHL 73 spectra to a three Gaussian model. Symbols denote the 
so-called BLVC (open dots), RVC (filled dots) and BHVC (filled triangles) 
velocity components (see text). 



with values increasing towards the south. At the peak of the emis- 
sion, the BHVC reaches a value of —100 km s _1 . The velocities 
of the BLVC range from -50 to -2 km s -1 , with -20 km s" 1 at 
the peak of the emission, with values increasing towards the east. 
Finally, the velocities of the RVC range from +20 to +40 km s _1 , 
with +35 km s" 1 at the peak of the emission. The highest values 
of the RVC are found around the emission peak, and the velocity 
decreases towards the edges of HHL 73. 

The line-profile decomposition also shows significant differ- 
ences, both in the spatial distribution and in the relative strength of 
the emission associated with each of the three velocity components 
(see Fig. [16] where the flux contours of each kinematic component 
are overlaid on the [S Il]-band NOT image). Regarding the spatial 
distribution, the BLVC and RVC components are present in almost 
all the spectra, while the BHVC component is concentrated in a 
smaller region, west of the emission peak. Regarding the relative 
strength of the components, the emission from the BHVC is the 
weakest (with flux values roughly one third of those found in the 
other components). The strongest emission comes from the BLVC 
component. Finally, the emission with the widest spatial distribu- 
tion comes from the RVC, the strength of its emission being slightly 
weaker than that of the BLVC. 

Interestingly, we found a significant spatial displacement 
among the positions of the emission peaks of the three velocity 
components (see also Figure |T6t. The emission peaks of the BLVC 
and RVC components are closely coincident (the estimated offset 
between the peaks of both components is ~ 1.3 arcsec; note that the 
IFS pixel size is 2.16 arcsec). In contrast, the BHVC peaks ~ 3.7 
arcsec northwest of them. As can be seen in Fig. [16] the emission 
from the redshifted component closely reproduces the spatial distri- 
bution of the bright bow-shaped HHL 73 emission, while its contri- 
bution to the collimated, low-brightness emission tail is lower than 
those from the blueshifted components. The two blueshifted com- 
ponents appear mainly associated with the northwestern component 
of the bright bow-shaped HHL 73 emission and with the HHL 73 
collimated emission tail. The shape of the the BLVC emission ap- 
pears elongated in the east-west direction, extending towards the 
west, and its emission peak closely coincides with the peak of the 
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Figure 15. Observed spectra (black histograms) and fits (dashed magenta 
lines) obtained from the model of three Gaussian components applied to 
the [S II] AA6717, 6731 A emission lines. The individual components are 
shown at the bottom of each panel. The same colour indicates the same 
component in each line. Residuals obtained subtracting the fit from the ob- 
served spectra are plotted in green. Labels in the left corner of each plot 
indicate the offset of each observed spectrum from the [S II] emission peak. 



[S II] emission found in the narrow-band NOT images. Finally, the 
emission from the BHVC, located in a smaller area, also appears 
slightly elongated in the east-west direction, peaking ~ 4 arcsec 
west of the [S II] emission peak. Note however that a clear asso- 
ciation of each kinematic component with a given morphological 
component of the HHL 73 [S II] emission cannot be easily de- 
rived only by looking to the [S II] channel maps, since the emis- 
sion from the wings of the stronger BLVC and RVC emissions give 
an appreciable contribution to the emission at positions where the 
BHVC (the kinematics component with the weakest emission) is 
located (e. g. west of the peak of the emission), thus contaminating 
the global emission at these positions. Finally, it should be pointed 
out that similar kinematic structures consisting of two blueshifted 
(low and high ) velocity componen ts have already been detected in 
Class I YSO |Takami et alj|2006j) and in microjets from T Tauri 
stars JBacciotti et alll200Cl : lLavallev-Fouquet. Cabrit & Dougadosl 
|2002b . 



4 DISCUSSION 

4.1 The origin of the HHL 73 optical emission 

In all the positions covered by the HHL 73 emission, the spectra 
extracted from the IFS datacube appear as low-ionization, line- 
emission spectra with no appreciable continuum emission, (see 
Figsll2landll3t. In particular, the spectrum shown in Fig.Qj] corre- 
sponding to the emission of the ~ 2 arcsec region around the peak 
position, shows no significant continuum emission superposed on 
the observed lines. 

We compared the [N II]/[S II] and [N Il]/Ha line ratio values 




Figure 16. Top-left: [S II] intensity maps obtained from the fit of the three- 
Gaussian model to the emission line profiles and adding the individual 
fluxes obtained for each velocity component. Top-right: [S II] intensity map 
for the BLVC. Bottom-left: [S II] intensity map for the RVC. Bottom-right: 
[S II] intensity map for the BHVC. The intensity of each velocity compo- 
nent is overlaid on the narrow-band [S II] NOT image. 



measured at several positions around the peak emission (where the 
signal-to-noise ratio is high enough to do the diagnostic with con- 
fidence) with those obtained from the diagnostic diagrams derived 
from the models of iKewlev & Dopital ( 120021) . These models give 
the ionization parameter (i. e. ionizing photon flux through an unit 
area) of the photoionized gas found in Hll and star-forming re- 
gions as a function of the line-flux ratios for a set of metallicities. 
The values found for the HHL 73 [N n]/[S II] and [N II]/Hq line ra- 
tios cannot be simultaneously fitted by the photoionization models: 
the [N Il]/Ha ratios are compatible with a low ionization parameter 
for metallicities z ^ 0.5 zq\ in contrast, for a wide range of ion- 
ization parameters, the [N II]/[S II] ratios are only compatible with 
metallicity ranges significantly lower (0.2 ^ z ^ 0.5 Zq). Thus, it 
does not seem plausible to interpret the HHL 73 optical emission 
as mainly originating from photoionization from a stellar source. 

In contrast, the [N II]/Hq and [S II]/Hq line ratio values de- 
rived in HHL 73 are fully compatible with those found for HH 
objects (i. e. from shock-excite d emission) i n the current diagnos- 
tic diagrams such as those of ICantdl jl98ll) . Both HHL 73 line 
ratio values correspond to diagram regions where the HH ob- 
jects are placed, being thus outside the regions where the spectra 
from photoionized nebula, such as the planetary nebula and the 
Hll regions, are located. It should be noted that the values mea- 
sured for these two line ratios (0.4—0.6 for [Nh]/Hq and 0.9-1 
for [S II]/Hq) are more compatible with the values found in HHs 
of high/intermediate excitation spectra than with the values found 
in HHs of low-excitation spectra. The degree of excitation of the 
HHL 73 spectra cannot be properly settled without data from the 
emission in a bluer wavelength range including high excitation 
emission lines (e. g. [O III]). Note, however, that the HH objects in 
which centimetre continuum emission was reported currently show 
a high-excitation spectrum. 

Thus, the HHL 73 IFS data suggest that the emission arises 
mainly from shock-excited gas, characteristic of an HH object. 

The contribution from continuum emission, if any, should be 
very weak, as derived from our optical data. We generated contin- 
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uum images from the IFS datacube by integrating, at each position, 
the signal over the wavelength range free of lines emission. In these 
images, we have only been able to find hints of an emission en- 
hancement relative to the neighbouring image background around 
the HHL 73 position, but a at very low-level (< 3a). In addition, 
the narrow-band, off-line filter NOT image fail to detect continuum 
emission, although it has less depth than the IFS data. This putative 
continuum contribution could arise from light scattered by dust bet- 
ter than from a point-like stellar source, as it extends over several 
pixels. 

4.2 The origin of the near-//f emission associated with 
HHL 73 

Let us discuss the nature of the 2MASS source found near the peak 
position of HHL 73. 

Regarding its size, the 2MASS J-, H- and A"-band images show 
that the 2MASS source appears slightly extended when compared 
with other sources in the fields. Typical FWHM of the field point 
sources are 3.3 ± 0.2 arcsec, while the FWHMs of the 2MASS 
source are 5.7, 5.3 and 4.8 arcsec in the H- and K-bands respec- 
tively, corresponding to a deconvolved size of ~ 4 arcsec (~ 3500 
AU for a distance of 900 pc). 

Regarding the near-//? colour, the (J — H,H — K) colour- 
colour diagrams are commonly used to discriminate between pre- 
main seq uence objects and field (main sequence and giant) stars 
(see e.g. biha et al.1 [2004). The near-//? colour indices of the 
2MASS source derived from the 2MASS photometry, (H — K) = 
1.26 ± 0.09, (J - H) = 1.47 ± 0.11, place the source inside 
the area corresponding to YSOs of Class II, whose near-//? emis- 
sion originates both from a photosphere and a circumstellar disc, 
but its colours are also compatible, within the photometric uncer- 
tainties, with being a Class I source, having thermal emission from 
a circumstellar envelope. Thus, from the near-//? colours, it seems 
plausible to infer that the 2MASS source corresponds to a very red, 
embedded Class II or a Class I source (i. e. a typical outflow-driving 
source ). According to current evolutionary models of PMS stars 
(e. g. IPalla & Stahleil 1 19931 ; ISiess, Dufour & ForestinJ |2000|) . the 
2MASS source lies in the region of the (J, J—H) colour-magnitude 
diagram that corresponds to low-mass YSOs with masses ^ LMq. 
The size derived for the 2MASS source is compatible with that ex- 
pected for an evacuated cavity around a YSO. However, in this sce- 
nario the bulk of the near-//? emission should mainly arise from 
continuum scattered sta rlight, and this does not seem to be the case. 
Davi s et al. I d200lLl2003b obtained high-resolution long-slit spectra 
at the HHL 73 position and were unable to detect continuum emis- 
sion in the H and K bands including [Fell] (at 1.64 pirn) and H2 
(at 2.12 /im) emissions. Thus, the continuum near-//? emission in 
HHL 73, if present, is weak. 

An alternative origin for the bulk of the near-//? emission 
is shock emission instead of scattered starlight. In this case, the 
2MASS Zf-band luminosity would mainly arise from the H2 2.12 
/im line, the //-band luminosity from the [Fell] 1.64 /im line, and 
the ./-band luminosity from other lines, like the [Fell] 1.25 /im 
line. However, this is hard to confirm, since no spectrum covering 
the ./-band range is available, and the 2MASS K and H magnitudes 
cannot be compare d reliably with the magnit udes derived from the 
line fluxes given bv lDavis et alj yOOl, 2003) - the line fluxes cor- 
respond to the integrated emission within a 0.8 x 2.7 arcsec, i. e. 
only a small fraction of the source size. In addition, the colour in- 
dices derived for the 2MASS source are also compatible wit h those 
expected for reddened molecular shocks, as was modelled by Smith 



dl995h for a wide range of shock parameters of molecular C- and 
J-shocks, who found that embebded shocks could mimic the near- 
//? colours derived for Class sources, both sh aring the same area 
of the near-//? colour diagrams. Furthermore, iDavis et all teOOll 
2003), from the analysis of the kinematics derived from the near- 
//? spectra, suggest that the [Fe II] and H2 emissions are tracing an 
outflow lying close to the plane of the sky, with a blueshifted com- 
ponent towards the west and a redshifted component towards the 
east. 

Thus, we found no conclusive evidence for the bulk of the 
near-//? emission of the 2MASS source close to the HHL 73 peak 
to be mainly tracing the scattered light from the YSO driving the 
outflow. Alternatively, a shock origin for the bulk of the near-//? 
emission seems to be compatible with the observations. 



4.3 The origin of the centimetre emission associated with 
HHL 73 

Centimetre continuum emission arising fro m HH objects has been 
detect ed in both high-mass (e. g. HH 80-8 1: lRodriguez & Reipurfr] 
19891: [M arti. Ro driguez & Reipurfhll 19931. a nd IRAS 16547-4247: 
Garav et al.1 120031 ; iRodriguez et al.1 120051) and low-mas s star- 
forming regions (e. g . HH 1 and HH 2: IPraydo et al.1 1 19851 : 
iRodrfguez et aill990l and HH 32A: lAnglada et alJll992l) Tln high- 
mass star-forming regions, the centimetre emission has a negative 
spectral index, which is consistent with non-thermal synchrotron 
emission arising from shocks resulting from the interaction of a 
collimated stellar wind with the surrounding medium, while in 
low-mass star-forming regions the centimetre emission shows a flat 
spectrum, which is consistent with optically thin free-free thermal 
emission from ionized gas. 

The spectral index of the centimetre HHL 73 emission can- 
not be derived from the present single-frequency observations. Fur- 
ther observations at another frequency are needed to estimate the 
spectral index in order to discriminate between a thermal and non- 
thermal origin for the centimetre emission of HHL 73. However, 
the exciting source of HHL 73 should be a low-/intermediate-mass 
YSO, since its bolometric luminosity is low (see next section). 
Thus, the centimetre emission detected in HHL 73 is most prob- 
ably free-free thermal emission from ionized gas. 

Two mechanisms could ionize the gas: photoionization or 
shocks. The observed flux density (Table [TJ, at the distance of 
HHL 73, corresponds to a c entimetre continu um luminosity of 0.21 
mjy kpc 2 . As discussed by I Angladal dl995h . the bolometric lumi- 
nosity required to produce such a centimetre continuum luminosity 
by photoionization is ~ 10 3 Lq, more than an order of magni- 
tude higher than the upper limit derived for the HHL 73 exciting 
source. In contrast, the HHL 73 centimetre continuum luminosity 
is typical of thermal radio-jets tracing molecular-outflow exciting 
sources, whose radio-continuum emission is produced by shock 
ionization. In addition, the momentum r ate of the HHL 73 m olecu- 
lar outflow, derived from the CO data of iDobashi et al.1 H993I) . fits 
the momentum rate vs. centimetre continuum luminosity relation- 
ship expected for shock-ionization induced e mission in mole cular- 
outflow sources of low bolometric luminosity (Anglada 199^). 

In conclusion, the centimetre continuum emission of HHL 73 
most probably arises from shock-ionized gas. 
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4.4 Is IRAS 21432+4719 tracing the exciting source of 
HHL 73? 

IRAS sources in star-forming regions are frequently found tracing 
the driving source of molecular and atomic outflows. The IRAS 
source closest to HHL 73 is 21432+4719, whose position from 
the IRAS Point Source Catalogue is ~ 30 arcsec northeast of the 
HHL 73 peak position (~ 0.13 pc for a distance of 900 pc), with a 
large uncertainty, so that HHL 73 lies outside the IRAS error ellipse, 
near its southwestern border (see Fig.[TJi. 

IRAS 21432+4719 is well detected in three IRAS bands, hav- 
ing a rising far-//? spectral energy distribution, with flux densities 
of 1.42, 12.19 and 23.45 Jy at 25, 60 and 100 pm respectively, 
and an upper limit of 0.25 Jy at 12 pm. The bolometric luminos- 
ity derived from the I RAS flux densities, following the recipe of 
IConnellev etal] J2007I) . is 27 Lq. This luminosity corresponds to 
a ZAMS star of 2.6 Mq. The IRAS colours derived from the fluxes 
are typical of embed ded protostellar sources l iBeichman et alj 19861 : 
iDobashi et alll 1992b. being redder than those found for T Tau-like 
sources ( Harri s, Clegg & Hugheslll988l) . 

Let us discuss whether IRAS 21432+4719 is tracing the em- 
bedded exciting source of HHL 73. 

As can be seen in Fig.Q] HHL 73 lies close to the border of a 
region of high visual extinction. A detailed study of the brightness 
distribution of the optical and near-//? emission of HHL 73 shows 
that there is a shift towards the northeast of the emission centroid 
position with increasing wavelength. The shift between the Ha and 
A"-band emission centroids is ~ 1.5 arcsec, while between the J- 
and H-, and between the H- and A"-band emission centroids it is ^ 1 
arcsec. These shifts are indicative of a large gradient of extinction, 
with extinction increasing towards the northeast. Thus, if there is 
far-//? emission associated with HHL 73, we expect to find it offset 
to the northeast of the optical and near-//? emissions. 

As already mentioned, the catalogue position of 
IRAS 21432+4719 is offset from HHL 73 in the expected di- 
rection, although somewhat displaced from the HHL 73 position. 
Nonetheless, the region encompassed by the IRAS error ellipse, i. e. 
the region with the highest probability of finding the IRAS source, 
almost reaches the HHL 73 location. Thus, the far-//? emission 
could be located nearer to the optical and near-//? emission than 
the IRAS catalogue position. 

No source is detected in the area encompassed by the IRAS 
error ellipse in the four bands of the Midcourse Space Experiment 
(MSX), at wavelengths from 8 to 22 pm. However, it should be 
noted that at 8.28 ^m, in the highest-sensitivity MSX band, there 
are hints of an enhancement in emission over the local background 
in a few pixels covering the position of the 2MASS source, but at a 
level below 3 times the image rms. 

In conclusion, in spite of the angular separation between the 
IRAS source and the optical and near-//? emission of HHL 73, 
IRAS 21432+4719 is the most plausible candidate to be driving 
HHL 73, although the evidence is inconclusive and new observa- 
tions in the millimetre wavelength range are needed to confirm this 
issue. 



4.5 A speculative scenario for HHL 73 

The HHL 73 m orphology is reminiscent of the HH 83 system (e. g. 
lReipurth|[l98^) . IRAS 05311-0631, the proposed exciting source 
of HH 83, is found located within a region with high visual ex- 
tinction, close to the apex of the optical emission. The HH 83 sys- 
tem is interpreted as an optical jet emerging from a conical cav- 
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Figure 17. Sketch of the scenario proposed for HHL 73. 



ity evacuated by the outflow. The counterjet is not visible since it 
flows inside the high-extinction region. The bulk of the bow-shaped 
emission at the base of the HH 83 jet arises from reflected light of 
the driving source, tracing the walls of the cavity. Apparently, the 
HH 83 system could mimic the HHL 73 scenario. 

However there are significant differences between both ob- 
jects. First, the bulk of the HHL 73 emission detected is line emis- 
sion, with small, if any, continuum contribution. In contrast, the 
emission from the bow-shaped base of the HH 83 jet has an appre- 
ciable continuum contribution. Second, the radial velocities mea- 
sured in HHL 73 present both redshifted and blueshifted values, 
while the radial velocities along the HH 83 jet are only blueshifted. 
Thus, a different scenario with another geometry than that proposed 
for the HH 83 system should be assumed to explain the spatial dis- 
tribution of the HHL 73 velocities. In addition, the separation be- 
tween IRAS 21432+4719 and the optical and near-//? emission of 
HHL 73 is appreciably larger than that between IRAS 0531 1-0631 
and the base of the HH 83 jet (~ 13 arcsec or ~ 0.03 pc for a dis- 
tance of 460 pc). However, as already discussed, IRAS 21432+4719 
could be located closer to the apex of HHL 73 than the position 
given by the IRAS Point Source Catalogue. 

Let us describe a qualitative scenario compatible with the 
complex HHL 73 kinematics derived from IFS data. Let us assume 
that the YSO exciting HHL 73 is deeply embedded and located to 
the northeast of HHL 73, near the apex of the HHL 73 emission. 
This source is ejecting supersonic gas with a low-collimation pat- 
tern (as, e. g., from a disc- or X-wind) and, i n addition, is driving 
a highly-collimated outflow (see, for example. iTakami et alj (2006) 
and references therein, where several models giving rise to out- 
flows having low- and high- velocity components are discussed). 
The low-collimated gas flowing outside the dark cloud can give, by 
projection effects, both the low-blueshifted and the redshifted ob- 
served velocities (i. e. the main contributors to the BLVC and RVC), 
while the BHVC would mainly arise from the high-collimated out- 
flow. The global kinematic pattern obtained from such a scenario 
would then mimic that produced by an outflow having a redshifted 
and two (low and high) blueshifted velocity components, even if all 
the supersonic gas we are observing is actually flowing outside the 
dark cloud (see the sketch in Fig.ll7t, 
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5 CONCLUSIONS 

From the analysis carried out of the CCD and IFS HHL 73 obser- 
vations, complemented with those available from the literature we 
found: 

The bulk of the HHL 73 optical red emission is mainly line emis- 
sion from shock-excited gas, such as is usually found in HH ob- 
jects. Continuum emission from dust scattered starlight, if present, 
seems to be weak. 

High-spatial resolution, narrow-band images of HHL 73 re- 
solved its bow-shaped morphology into two components and a 
collimated emission ~ 4 arcsec long emerging towards the west, 
which is reminiscent of a microjet. 

The kinematics of HHL 73 derived from the IFS data show that 
the emission is supersonic, with both blueshifted and redshifted 
velocities. From a line-profile multicomponent model, we were 
able to resolve three velocity components for the [S II] AA6717, 
6731 A emission lines. These kinematic components correlate with 
the morphological components discovered from the high-resolution 
CCD images. The redshifted component, being the most extended 
one, spatially correlates with the bright, extended HHL 73 emis- 
sion, elongated in the north-south direction. Its spatial distribution 
is also in good agreement with that of the radio continuum emis- 
sion. The two blueshifted (low- and high- velocity) components are 
morphologically well correlated with the northwestern HHL 73 
emission and the microjet. In order to mimic the complex kine- 
matic pattern derived from the IFS spectra, we tentatively propose 
a scenario involving both low- and highly-collimated supersonic 
gas ejected from a YSO. 

In the light of the data here discussed, we propose that the ex- 
citing source of HHL 73 (and possibly of HH 379) is deeply em- 
bedded in the dark cloud causing the high extinction, lying to the 
northeast of HHL 73, close to its apex. In spite of its angular sep- 
aration, IRAS 21432+4719 is the most plausible candidate to be 
tracing the low-/intermediate YSO that is driving the centimetre, IR 
and optical emission associated with HHL 73, although high sen- 
sitivity observations with high angular resolution in the millimetre 
wavelength range, not affected by extinction, are needed to confirm 
this issue and to establish clearly the location and properties of the 
HHL 73 driving source. 

Further observations, both in the centimetre wavelength range 
and IFS data covering the wavelength range blueward of Ha, to 
include high-excitation HH emission lines, are needed to settle the 
nature of the HHL 73 object completely. 
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